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Abstract 

We consider an undirected graph G — {VG, EG) with a set T C VG of terminals, 
and with nonnegative integer capacities c{v) and costs a[v) of nodes v € VG. A 
path in G is a T-path if its ends are distinct terminals. By a multiflow we mean 
a function F assigning to each T-path P a nonnegative rational weight F{P), and 
a multiflow is called feasible if the sum of weights of T-paths through each node v 
does not exceed c{v). The value of F is the sum of weights F{P), and the cost of F 
is the sum of F{P) times the cost of P w.r.t. a, over all T-paths P. 

Generalizing known results on edge-capacitated multiflows, we show that the 
problem of finding a minimum cost multiflow among the feasible multiflows of maxi- 
mum possible value admits half-integer optimal primal and dual solutions. Moreover, 
we devise a strongly polynomial algorithm for finding such optimal solutions. 
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1 Introduction 



1.1 Multiflows 

For a function (p: X ^ M4, and a subset A C X, we write '^{A) to denote ^^eA'-Pi-'^)- 
The incidence vector of A in R'''" is denoted by x^-, i-^- X^(c) is 1 for e A and for 
e & X — A (usually X will be clear from the context). When j4 is a multiset, x^i^) 
denotes the number of occurrences of e in yl. 

In an undirected graph G, the sets of nodes and edges are denoted by VG and EG, 
respectively. When G is a directed graph, we speak of arcs rather than edges and write 
AG instead of EG. A similar notation is used for paths, cycles, and etc. 

A walk in G is meant to be a sequence {vq, ei,vi, . . . , ek,Vk), where each Cj is an edge 
(or arc) and Vi-i,Vi are its endnodes; when G is a digraph, e, is directed from Vi-i to Vi. 
Edge-simple (or arc-simple) walks are called paths. 

We consider an undirected graph G and a distinguished subset T C VG of nodes, 
called terminals. Nodes in VG — T are called inner. A T-path is a path P in G whose 
endnodes are distinct terminals; we usually assume that all the other nodes of P are 
inner. The set of T-paths is denoted by V. A multiflow is a function F: V Q-|-. 
Equivalently, one may think of F collection 

(1.1) {(ai,Pi),...,(a„,P„)} 

(for some m), where the Pi are T-paths and the are nonnegative rationals, called 
weights of paths. Sometimes (e.g., in |IKN98| ) such a multiflow F is called free to 
emphasize that all pairs of distinct terminals are allowed to be connected by flows. The 
value val(-F) of F is F{P). For a node v, define 

(1.2) F(v) := ^(F(P): P G E l^P); 

the function F on VG is called the ( node ) load function. Let c : VG — )• be a nonneg- 
ative integer function of node capacities. We say that F is feasible if F{v) < c{v) for all 
V G VG. 

Suppose we are given, in addition, a function a: VG — > Z+ of node costs. Then the 
cost a{F) of a multiflow F is the sum a{P)F{P), where a{P) stands for the cost 
a{VP) of a path P. 

In this paper we consider the following problem: 

(N) Given G, T, c, a as above, find a multiflow F of minimum possible cost a{F) among 
all feasible multiflows of maximum possible value. 

1.2 Previous results 

When |r| = 2, (N) turns into the undirected min-cost max-flow problem under node 
capacities and costs, having a variety of applications; see, e.g., |FF62l lLa76| . It admits 
integer optimal primal and dual solutions [FF62J. 

In the special case a = 0, we are looking simply for a feasible multiflow of maximum 
value. Such a problem has half-integer optimal primal and dual solutions, due to results 
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of Pap |Pa07| and Vazirani (VaOl) . respectively. Also it is shown in [Pa07j that the 
problem is solvable in strongly polynomial time by using the ellipsoid method. 

An edge-capacitated version of (A^) has been well studied. In this version, denoted 
by (E), c and a are functions on EG rather than VG. For a multiflow F, its edge load 
function is defined similarly to (|1.2p : 

(1.3) F{e) := ^ (F(P) : Per,ee EP) for all e G EG, 

and its cost is defined to be a{EP)F{P). Problem (E) is reduced to (A^) by adding 
an auxiliary node on each edge, but no converse reduction is known. 

An old result is that (E) has a half-integer optimal solution |Ka79| . Also it is shown 
in |Ka94| that (E) has a half-integer optimal dual solution and that half-integer primal 
and dual optimal solutions can be found in strongly polynomial time by using the ellip- 
soid method. A "purely combinatorial" weakly polynomial algorithm, based on cost and 
capacity scaling, is devised in |GK97| . 

In the special case of (E) with a = 0, the half-integrality results are due to 
Lovasz |Lo76| and Cherkassky |Ch77) . and a strongly polynomial combinatorial algo- 
rithm is given in |Ch77| (see also [IKN98| for faster algorithms). 

1.3 New results 

In this paper we prove that (N) always admits a half-integer optimal primal and dual 
solutions. In particular, this implies all half-integrality results mentioned in the previous 
subsection. 

Similar to |Ka94| . we introduce a parametric generalization of (A^), study properties 
of geodesies (shortest T-paths with respect to some length function), and reduce the 
parametric problem to a certain single-commodity flow problem. However, the details of 
this construction are more involved. In particular, the reduced problem concerns integer 
flows in a hidirected graph. 

The second goal is to explore the complexity of (A^). We show that half-integer 
optimal primal and dual solutions to the parametric problem (and therefore to (A^)) can 
be found in strongly polynomial time by using the ellipsoid method. 

2 Preliminaries 

2.1 Parametric problem and its dual 

Instead of (A^), it is convenient to consider a more general problem, namely: 

(A^a) Given G,T,c,a as in (N ) and, in addition, A € find a feasible multiflow F 
maximizing the objective function ^{F, a. A) := A • val(-F) — a{F). 

We will prove the following 

Theorem 2.1 For any A € "Lj^, problem (N\) has a half-integer optimal solution. 

(Note that ^{F, qa, qX) = q ■ ^{F, a, A) for any multiflow F and q £ Q+. Therefore, the 
optimality of a multiflow in the parametric problem preserves when both a and A are 
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multiplied by the same positive factor q. This implies that the theorem is generalized 
to arbitrary a : VG Q+ and A € Q-|- (but keeping the integrality of c). However, we 
prefer to deal with integer- valued a and A in what follows.) 

By standard linear programming arguments, {N) and (A^a) become equivalent when 
A is large enough (moreover, the existence of a half-integer optimal solution for (A^a) 
easily implies that taking A := 2ciyG)a{yG) -|- 1 is sufficient). 

Problem {N\) can be viewed as a linear program with variables F{P) G Q^- assigned 
to T-paths P. Assign to a node v € VG a variable l[v) G Q-|-. Then the linear program 
dual to (A^a) is: 

(Da) Minimize c ■ I provided that the following holds for every T-path P: 
(2.1) l{P)>\-a[P). 

2.2 Translating to edge lengths 

The above dual problem {D\) involves lengths of paths (namely, 1{P) and a{P)) deter- 
mined by "lengths" of nodes {I and a, respectively). It is useful to transform lengths of 
nodes into lengths of edges. To do so, for w: VG — t- Q^., we define the function w on 
EG by 

(2.2) W{e) := auw{u) + a^wlv) for e = uv € EG, 

where := ^ if a; € VG — T, and := 1 if x € T. This provides the correspondence 

(2.3) w{P) = w{P) for each T-path P 

(where w{P) stands for w{VP), and w{P) for w{EP)). For a, / as above, define 

(2.4) £:=l + a. 

Let dist£(u, f) denote the ^-distance between vertices u and f, i.e. the minimum -^-length 
^{P) of a u-v path P in G. Then, in view of ()2.3p and (12. 4p . the constraints in {D\) can 
be rewritten as 

(2.5) dist£(s, t) > A for all s,teT,Sy^t. 

By the linear programming duality theorem applied to {Nx) and (Dx), a feasible 
multiflow F and a function /: VG Q+ satisfying (j2.5p are optimal solutions to (A^a) 
and (Dx), respectively, if and only if the following (complementary slackness conditions) 
hold: 

(2.6) if P is a T-path and F{P) > 0, then i{P) = A; in particular, P is ^-shortest; 

(2.7) if w S VG and l{v) > 0, then v is saturated by F, i.e. F{v) = c{v). 

In the rest of the paper, to simplify technical details, we will always assume that the 
input costs a of all nodes are strictly positive. Then the edge lengths i defined by ()2.4p 
are strictly positive as well. This assumption will not lead to loss of generality in essence, 
since the desired results for a nonnegative input cost function a can be obtained by 
applying a perturbation technique in spirit of |Ka94t pp. 320-321] (by replacing a by an 
appropriate strictly positive cost function). 
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2.3 Geodesies 



Condition (|'2.6|) motivates the study of the structure of ^-shortest T-paths in G. To this 
aim, set p := min{dist£(s, t) | s, t G T, s 7^ t}. A T-path P such that i{P) = p is caUed 
an i-geodesic (or just geodesic if i is clear form the context). When a multiflow F in G 
is given as a collection (jl.ip in which all paths Pi are ^-geodesies, we say that F is an 
i-geodesic multiflow. 

Next we utilize one construction from |Ka79| [Ka94j . with minor changes. Consider 
a node v G VG. Define the potential it{v) to be the ^-distance from v to the nearest 
terminal, i.e. tt{v) := minjdist^ (f , t) | t G T}. Set VGg := {v G VG \ tt{v) < ^p} (in 
particular, T C VGi). For s G T, define := {v G VG \ dist^(s,r;) < ^p}- Also define 

:= {v eVG\ tt{v) = ip}. We refer to V as the zone of a terminal s € T, and to V'^ 
as the set of central nodes (w.r.t. i). The sets (s G T) and V'^ are pairwise disjoint 
and give a partition of VG^. 

The following subset of edges is of importance: 

EGe := {uv e EG\3s £T : ueV\ve V UV\\7r{u) - 7r{v)\ = £{uv)} 
U{uv e EG \3s,t eT,s ^t : u eV'',v eV\ Tr{u) + -k{v) + e{uv) = p}. 

One can see that the subgraph Gi := {VGe, EGe) of G contains all ^-geodesies. 
Moreover, a straightforward examination shows that the structure of geodesies possesses 
the properties as expressed in the following lemma (which is, in fact, a summary of 
Claims 1-3 from jKa94j and uses the strict positivity of £). 

Lemma 2.2 Let P be an i-geodesic running from s G T to t (zT. Then P is contained 
in Gi and exactly one of the following takes place: 

1. P contains no central nodes and can he represented as the concatenation Pi o 
(ti, e, v) o P2, where u G , v G V^, s ^t, and e G EG^. 

2. P contains exactly one central node w ^V'^ and can he represented as the concate- 
nation P = Pi o (li, ei,w, e2,v)oP2, where u G V^ , v G V^, s ^ t, and ei, 62 G EGe. 

In hoth cases, parts Pi and P2 are contained in the induced suhgraphs Ge[V^] and Ge[V*], 
respectively. The potentials it are strictly increasing as we traverse Pi from s to u, and 
strictly decreasing as we traverse P2 from v to t. 

Conversely, any T-path in Ge oheying the ahove properties is an i-geodesic. 

3 Primal half-integrality 
3.1 Auxiliary bidirected graph 

In this subsection we introduce an auxiliary bidirected graph, which will be the cor- 
nerstone of our approach both for proving half-integrality results and for providing a 
polynomial-time algorithm. 

Given G, T, c, a and A as above, let I be an optimal solution to (Dx). Form the edge 
lengths i :=1 -\-a, the potential tt, the subgraph Ge, and the sets {s G T) and V\ as 
in Subsection 12.31 One may assume that p := min{dist£(s, i) | s, t G T, s 7^ t} = A (since 
J3 > A, by (|2.5p . and if p > X then P = is an optimal solution to {N\), by ()2.6p ). 
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For further needs, we reset c := 2c, making all node capacities even integers. Now our 
goal is to prove the existence of an integer optimal multiflow F in problem {N\) (which 
is equivalent to proving the half-integrality w.r.t. the initial c). 

Recall that in a bidirected graph (or a BD-graph for short) edges of three types 
are allowed: a usual directed edge, or an arc, that leaves one node and enters another 
one; an edge directed from both of its ends; and an edge directed to both of its ends 
(cf. |EJ70l[Sc03| ). When both ends of an edge coincide, the edge becomes a loop. For our 
purposes we admit no loop entering and leaving its end node simultaneously. Sometimes, 
to specify the direction of an edge e = uv at one or both of its ends, we will draw arrows 
above the corresponding node characters. For example, we may write ii^ if e is directed 
from u to V {a usual arc), 17"V if e leaves both u, v, if e enters both u, v, and itv if 
e leaves u (and either leaves or enters v). 

A walk in a BD-graph is an alternating sequence 

P = {s = vo,ei,vi, . . . ,ek,Vk = t) 

of nodes and edges such that each edge Cj connects nodes and Vi, and for i = 

1, . . . ,k — 1, the edges Cj, Cj+i form a transit pair at Vi, which means that one of Cj, Cj+i 
enters and the other leaves fj. As before, an edge-simple walk is referred to as a path. 

Now we associate to Gi a BD-graph H with edge capacities c: EH — > Z_|_, as follows 
(see Fig. [1] for an illustration). Each noncentral node v € VGi — Y'^ generates two nodes 
in }i. They are connected by edge (arc) e^, going from to and having the 
capacity equal to c(f). We say that inherits the capacity of the node v. For s €z T, 
the set V := {v^,v'^ \ v G V^} in H is called the zone of s, similar to in G. 

Consider an edge e = uv EGi. Let u,v G for some s € T and assume for 
definiteness that 7r(n) < 7r(f) (note that £{uv) > implies 7r(n) ^ vr(u); this is where 
the strict positivity of the cost function a is important). Then e generates in H an edge 
(arc) going from to v^, and we assign infinite capacity to it. (By "infinite capacity" 
we mean a sufficiently large positive integer.) Now let u (z and v (z V* for distinct 
s,t Then e generates an infinite capacity edge it'^V^ (leaving both u'^ and v'^). 

The transformation of central nodes is less straightforward. Each w G generates 
in H a so-called gadget, denoted by F^. It consists of \T\ + 1 nodes; they correspond to 
w and the elements of T and are denoted as 9^ and 9w^s, s G T. The edges of F^ are: a 
loop leaving 9w (twice) and, for each s £ T, an edge e^^^ going from 9w,s to 9w, called 
the s-leg in the gadget. Each edge in F^ is endowed with the capacity equal to c{w). 

Each gadget F^ is connected to the remaining part of H as follows. For each edge 
of the form vw in Gi, we know that v & for some s (z T (by the construction of Gg). 
Then vw generates an infinite capacity edge (arc) going from v'^ to 9^ 

Finally, we add to H an extra node q, regarding it as the source, and for each s G T, 
draw an infinite capacity edge (arc) from q to s^. 

The obtained BD-graph H captures information about the ^-geodesies in G. Namely, 
each £- geodesic P going from s to t induces a unique closed q-q walk P \n H. The first 
and the last edges of P are and t^% , respectively. For a noncentral node v in P, P 
traverses the edge if^lf'^. An edge uv G EP with 7r(ii) < Tr{v) inside a zone induces the 
edge if^lt^ in P. An edge uv G EP connecting different zones (if any) induces the edge 
it^tr^ in P. Finally, suppose P traverses a central node w € and let uw, wv G EGg 
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(a) Graph Gt- (b) Bidirected graph H. 

Figure 1: Constructing graph H. Here T = {p, s,r}, = {p, a,6}, = {s,c, d}, 
y = {r, e, /,(7}, V'^ = {w}. (The source q is not shown.) Bidirected edges leaving one 
endpoint and entering the other are indicated by ordinary directed arcs. Marked are one 
f-geodesic P and its image P. 

be the edges of P incident to w. By Lemma 12.21 n € ^ and v £ for some s ^ t. 
Then the sequence of nodes u, w,v m P generates the subpath in P with the sequence of 

edges u O^^^g^ c^^^s^ ^Wjt"! ^w.t'^ • 

The resulting walk P is edge-simple, so it is a closed path. Conversely, let Q be 
a (nontrivial) q~q walk in H. One can see that Q with q removed is concatenated as 
Qi o Q' o Q2, where Qi is a directed path within a zone V , Q2 is reverse to a directed 
path within a zone V (with possibly s = t), and Q' either (i) is formed by an edge 
T^^tr^ connecting these zones (in which case s 7^ t), or (ii) is the walk with the sequence 
of edg 6S dw ,s ^ (^w : (^w ,t 1 for some central node w of Gi- Moreover, the image in G of each 
of Qi,Q2 is an ^-shortest path. When s = t happens in case (ii), Q traverses the edge 
Gw,s twice. In all other cases, Q is edge-simple and its image in G is an ^-shortest T-path 
(a A-geodesic). 

These observations show that there is a natural bijection between the i'-geodesics in 
G and the (nontrivial) q-q paths in H. 

We will refer to the BD-graph H described above as the compact BD-graph related to 
Gi] it will be essentially used to devise an efficient algorithm for solving (A^a) hi Section[H 
Besides, in the proof of the primal integrality (with c even) in Section [5l we will deal with 



7 



a modified BD-graph. It is obtained from H as above by replicating each gadget into 
c{w) copies r^i, i = 1, . . . ,c(t(;), cahed the 1-gadgets generated by w. More precisely, 
to construct T^i, we make i-th. copy of the node 6^, i-th copy e^i of the loop 
leaving 6*^1 (twice), and i-th copy e^i^^ of each leg ew,s, s € T, where e^i goes from 
6w,s to (so 6w,s, s (zT, are the common nodes of the created 1-gadgets). All edges in 
these 1-gadgets are endowed with unit capacities. 

We keep notation H for the constructed graph and call it the expensive BD-graph 
related to Gi. Also we keep notation c for the edge capacities in H. There is a natural 
relationship between the q-q walks (paths) in both versions of H. The 1-gadgets created 
from the same central node w of Gi are isomorphic, and for any i, j = 1, . . . , c{w), there 
is an automorphism of H which swaps and and is invariant on the other nodes. 

3.2 Bidirected flows 

Let r be a bidirected graph. Like in usual digraphs, 5™(u) and d"^^{v) denote the sets 
of edges in T entering and leaving v G VT, respectively. A loop e at v, if any, is counted 
twice in 5^^{v) if e enters v, and twice in 6°^^{v) if e leaves v; hence (5'°(t;) and S'^^^{v) 
are actually multisets. (Recall that we do not allow a loop which simultaneously enters 
and leaves a node.) 

Let g be a distinguished node with (5'"(g) = in F (the source) and let the edges of T 
have integer capacities c: ET — > A bidirected q-flow, or a BD-flow for short, is a 
function /: ET — ?> Q+ satisfying divj;(f) = for all nodes v € VT — {g}; and the value 
of / is defined to be div/(g) (cf. [ GK04| ). Here 

(3.1) div;(t;):=/(5°-*(t;))-/(5'"(^;)) 

is the divergence of / at v. Note that if e is a loop at v then e contributes ib2/(e) in 
d\Yf{v). If /(e) < c(e) for all e € ET then / is called feasible. In addition, if / is integer- 
valued on all edges then we refer to / as an integer bidirected g-flow, or an IBD-flow. 
One can see that finding a fractional (resp. integer) BD-flow of the maximum value is 
equivalent to constructing a maximum fractional (resp. integer) packing of closed q-q 
walks (they leave q twice). 

Return to an optimal solution / to {Dx), and let i := a + 1. Consider the (expensive 
or compact) BD-graph H related to Gi, and the capacity function c on the edges of 
H (constructed from the node capacities c of G). The above correspondence between 
^-geodesies in G and q-q paths in H is extended to ^-geodesic multiflows in G and 
certain g-flows in H (where q is the source in H as before). More precisely, let F be a 
(fractional) ^-geodesic multiflow in G represented by a collection of ^-geodesies Pi and 
weights := F{Pi), i = 1,. . . ,m (cf. (jl.ip ). Then each Pi determines a q-q path Pi 
in H, and / := aix^^^ + . . . + UmX^^"^ is a BD-flow in H; we say that / is generated 
by F (note that val(/) = 2val(F)). Furthermore, / is feasible if F is such, and for each 
central node w €V\ the following relations hold: 

(3.2) ^/(eu.,s) = 2/(e^) and /(e^,^) < /(e«,) for each s G T. 

Considering an arbitrary BD-flow / in H, we say that / is good if it satisfies (|3.2p for 
all w [h ere the second relation in (|3.2|) is important, while the first one obviously 

holds for any BD-flow). The following assertion is of use. 
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Lemma 3.1 Let f be a good BD-flow in H. Then f is generated by an (.-geodesic 
multiflow F in G. Moreover, if f is integral, then it is generated by an integer (.-geodesic 
multiflow F. In both cases, F can be found in 0{\EH\) time. 

Proof. Suppose there is a central node w ^ V''^ such that /(e^) > 0. Let us say that 
p dominates at w (w.r.t. /) if f{eyj^p) = /(e^,). From ()3.2p it follows that there exist 
distinct s,t ^ T such that /(e^^^j), f{ew,t) > and none of p G T — {s, t} dominates at w. 
Choose such s, t. Build in a maximal walk Q starting with 0^, eu,,s,dw,s, ■ ■ ■ and such 
that /(e) > for all edges e of Q. It is easily seen from the construction of H that Q 
is edge-simple, terminates at q, and have all vertices in V , except for Oyj^Oyj^g- Build a 
similar walk (path) Q' starting with Oyj,ewfi,Oyj^f Then the concatenation of the reverse 
to Q, the loop Cw and the path Q' is a q-q path and its image P in G is an geodesic 
(from s to t). 

Assign the weight of P to be the maximum number a subject to two conditions: 
(i) Oi < /(e) for each e G EP, and (ii) the flow f '■= f — olx^ is still good. If a is 
determined by (i), we have |supp(/')| < |supp(/)| (where supp((/9) := {x \ ip{x) ^ 0}), 
whereas if a is determined by (ii), there appears p G T dominating at w (w.r.t. /'). If 
f'i^w) > 0, repeat the procedure for /' and w, otherwise apply the procedure to /' and 
another w' (zV\ and so on. (Note that if, in the process of handling w, a current weight 
a is determined by (ii), then the weights of all subsequent paths through e^ are already 
determined by (i); this will provide the desired complexity.) 

Eventually, we come to a good flow / with /(e^) = for all w G This / is 
decomposed into a sum of flows along q-q paths in a straightforward way, in 0{\EH\) 
time (like for usual flows in digraphs). Taking together the images in G of the constructed 
weighted q-q paths, we obtain a required ^-geodesic multiflow F. The running time of 
the whole process is 0{\EH\), and if / is integral, then the weights a of all paths are 
integral as well. (Integrality of a current weight a subject to integrality of a current flow 
/ is obvious when a is determined by (i), and follows from the fact, implied by (|3.2p . 
that for any p & T, ^^_^p f{ew,s) — f{&w,p) is even, when a is determined by (ii).) □ 

Remark 3.2 In case of an expensive BD- graph H , any feasible IBD-flow f is good. 
Indeed, for any 1-gadget T^i in H, we have c(e^i) = 1 and, therefore, f{ew) G {0)1}- 
The second relation in \3.2^) is trivial when /(e^) = 0, and follows from the constraints 
f{ew,s) < c(e^,s) = 1 (s gT) when /(e^) = 1. 

Define the following subset of edges in H: 

(3.3) Eo := {e^\v£ VG, l{v) > 0}. 

For an optimal (possibly fractional) solution F to {Nx) and a node v G VG with l(v) > 0, 
we have F{v) = c{v) (by (|2.7p ): so the edge e^ of H corresponding to v is saturated by 
the BD-flow / generated by F, i.e. /(e^,) = c(e„). We call the edges in Eq locked. 

Thus, the graph H admits a (fractional) good feasible BD-flow saturating the locked 
edges. The following strengthening is crucial. 

Proposition 3.3 There exists a good feasible IBD-flow in H that saturates all locked 
edges. 
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A proof of this proposition involves an additional graph-theoretic machinery and will 
be given in Section [5j Assuming its validity, we immediately obtain Theorem 12.11 from 
Lemma 13.11 

4 Solving (A^a) in strongly polynomial time 

In this section we devise a strongly polynomial algorithm for solving the primal parametric 
problem (Nx). As before, we assume that a and A are integral and that the node capacities 
c are even, so our goal is to find an integer optimal multiflow. 

The algorithm starts with computing a (fractional) optimal dual solution I and con- 
structing the BD-graph H w.r.t. the length function i := a + l. Then it finds a good IBD- 
flow f in H saturating the locked edges (assuming validity of Proposition 13. 3p . Applying 
the efficient procedure as in the proof of Lemma 13.11 to decompose / into a collection of 
paths with integer weights, we will obtain an integer optimal solution to (Nx). 

To provide the desired complexity, we shall work with H given in the compact form 
(defined in Subsection 13. 1|) . The core of our method consists in finding the load function 
of some integer optimal multifiow F in G (without explicitly computing F itself). This 
function will just generate the desired IBD-fiow in H. We describe the stages of the 
algorithm in the subsections below. 

4.1 Constructing an optimal dual solution 

Problem (Dx) straightforwardly reduces to a "compact" linear program, as follows. Be- 
sides variables l{v) € Q+ {v G VG), assign a variable <^s{v) S Q to each terminal s (z T 
and node v (a sort of "distance" of v from s). Consider the following problem (where I 
and a are defined according to (|2.2p ): 

(4.1) Minimize c • I subject to the following constraints: 

ips{u) - ips{v) < a(e) +lie) 

fsiv) — ips{u) < a(e) -|- 1(e) for each e = E EG; 
(fsit) — ^s{s) > A for all s,t & T, s ^ t. 

Lemma 4.1 Programs (Dx) and J^.il ) are equivalent. 

Indeed, if (/,(/?) is a feasible solution to (|4.ip then, obviously, I is a feasible solution 
to {Dx)- Conversely, let / be a feasible solution to {Dx)- For v G VG and s (zT, define 
(Ps{v) '■= disti{s,v), where £ ■.= l + a. It is easy to check that {l,yp) is a feasible solution 
to (liTjl . 

The size of the constraint matrix in (|4.ip (written in binary notation) is polynomial in 
Therefore, {Dx) is solvable in strongly polynomial time by use of Tardos's version 
of the ellipsoid method. (This remains valid when a and A are nonnegative rational 
numbers.) 
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4.2 Computing the load function of an optimal multiflow 

The following fact is of importance. 

Lemma 4.2 One can find, in strongly polynomial time, a function g: VG Z_|_ such 
that g{v) = F{v) for all v € VG, where F is some integer optimal multiflow in (Nx). 

Proof. We explain that in order to construct the desired it suffices to compare 
two optimal objective values: one for the original (integer) costs a and the other for 
certain perturbed costs a^. These values are computed by solving the corresponding dual 
problems by the method described in the previous subsection. 

More precisely, let vi,...,Vn be the nodes of G. Let U := maxjc(fj) + 2, define 
e{vi) := l/C/*"*"^, i = 1, . . . ,n, and define the cost function on VG to be a + e. Then 
for any integer feasible multiflow F, we have 

< a, A) - as, X) = Y1 F'{vi)e{vi) < U''^ + U'^ + . . . + < L 

i 

This and the fact that <J?(F, a, A) is an integer (as F,a,X are integral) imply that if F 
is optimal for a^, then F is optimal for a as well. (An integer optimal multiflow for 
even capacities c exists by Theorem 12.11 ) Moreover, for such an optimal F, the number 
r := ^iF{vi)e(vi) is computed in strongly polynomial time, since it is equal to c-l — c-lg, 
where I and l^ are optimal dual solutions for a and a^, respectively. Here we use the LP 
duality equalities $(-F, a, A) = c-l and ^{F, a^, A) = c-Z^. Also the size of binary encoding 
of Oe is bounded by that of a times a polynomial in n, so the dual problem with is 
solved in strongly polynomial time w.r.t. the original data. 

Hence, we have rU^~^^ = Xli The number rf/""*"^ is an integer and, in 
view of F{vi) < c{vi) < U for each i, the n coefficients in its base U decomposition (the 
representation via degrees of U) are just F(vi), . . . , F{vn), thus giving g. □ 

Recall that together with a node load function each multiflow F also induces its edge 
counterpart (see (|1.3p ). Lemma 14.21 can be strengthened as follows. 

Lemma 4.3 One can find, in strongly polynomial time, a function g: VG U EG TL^ 
such that g{v) = F{v) for all v € VG and g{e) = F{e) for all e € EG, where F is some 
integer optimal multiflow F in (Nx). 

Proof. Split each edge e = uv oi G into two edges uxg, XeV in series and assign to each 
new node Xg the capacity c(xe) := min{c(u), c(f )} and the cost a{xe) '■= 0. Clearly this 
transformation does not affect the problem in essence. The node load function, which 
can be found in strongly polynomial time by Lemma 14. 2[ yields the desired node and 
edge load functions in the original graph G. □ 

4.3 Constructing an optimal primal solution 

Now we explain how to find, in strongly polynomial time, an integer optimal multiflow 
solving (A'^^) (for a graph G, even node capacites c, rational node costs a, and an integer 
parameter A) by using an optimal dual solution I and a function g as in Lemma 14.31 For 
this g, there exists an integer ^-geodesic multiflow F' in G satisfying F'[v) = g{v) for all 



11 



V € VG and F'{e) = g{e) for all e € EG, where (. :=a + 1. Our goal is to construct one 
of such multiflows explicitly. 

To do this, we consider the compact BD-graph H related to Gi (see Subsection 13. ip 
and put / to be the function on EH corresponding to g. More precisely, let E' be the 
subset of edges of H neither incident to q nor contained in the gadgets Tyj {w € V'^). By 
the construction of H, there is a natural bijection 7 of E' to the set {VGi -V^)U EGe. 
For each e G E' , we set /(e) := g{'y{e)). In their turn, the values of / on the edges of a 
gadget are assigned as follows: for the loop at 9^,, set /(e^) := g{w), and for each 
s £ T, set f{ew,s) '■= J^idi^) '■ ^ ^ Eg^w), where Eg^w is the set of edges in Gi connecting 

and w. Finally, for each s € T, we set f{qs) := h{s). 

Using the fact that the function g on VG U EG is determined by some integer optimal 
(^-geodesic) multiflow F' , it follows that the obtained function / on EH is integer- valued 
and has zero divergency at all nodes different from q. So / is an IBD-flow in H. Moreover, 
/ is generated by F' as above; in particular, / is good (i.e. satisfies (|3.2p ). By Lemma [3.11 
we can find, in strongly polynomial time, an integer ^-geodesic multiflow F generating /. 
Then F and F' have the coinciding node and edge load functions, and the optimality of 
F' implies that F is an integer optimal solution to (Nx) as well, as required. 

5 Proof of Proposition 13.3 

To complete the proof of Theorem 12.11 it remains to prove Proposition 13. 3| which claims 
the existence of an IBD-flow saturating the "locked" edges. We eliminate the lower 
capacity constraints (induced by the locked edges) by reducing the claim to the existence 
of an IBD-flow with a certain prescribed value. 

5.1 Maximum IBD-flows 

Let r be a bidirected graph with a distinguished source q and edge capacities c : FT — )> Z+ , 
as described in Subsection 13.21 The classic max-flow min-cut theorem states that the 
maximum flow value is equal to the minimum cut capacity. A bidirected version of 
this theorem involves a somewhat more complicated object, called an odd barrier. In this 
subsection we give its definition and state the crucial properties (in Theorems 15 . 1 1 [5^ 15. 3p 
that will be used in the upcoming proof of Proposition 13.31 These properties are nothing 
else than translations, to the language of bidirected graphs, of corresponding properties 
established for integer symmetric flows in skew-symmetric graphs, as we will explain in 
the Appendix. 

Let X C VT — {q}. The flip at (the nodes of) X modifies F as follows: for each node 

V £ X and each edge e incident to f , we reverse the direction of e at u (while preserving 
the directions of edges at nodes in VG — X). For example, if e = and u,v G X 
then e becomes and e = ui and u ^ X 3 v then e becomes l^'V. BD-graphs F 
and F' are called equivalent if one is obtained by a flip from the other. Note that flips do 
not affect bidirected walks in F in essence and do not change the maximum value of an 
IBD-flow in it. We will essentially use flips to simplify requirements in the definition of 
odd barriers below. 

Next, we employ a special notation to designate certain subsets of edges. For X,Y CI 
VT, let [X, Y] denote the set of edges with one endpoint in X and the other in Y. We 
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Figure 2: A bidirected odd barrier. Grayed edges correspond to odd capacity constraints 

(w.r.t. r'). 



will often add arrows above X and/or Y to indicate the subset of edges in [X, Y] directed 
in one or another way. For example, [X , 1^] denotes the set of edges that enter both X 
and Y, , Y] denotes the set of edges leaving X and having the other endpoint in Y 
(where the direction is arbitrary), and ^Y^] denotes the set of edges that leave X at 
both endpoints (including twice leaving loops). When Y = VT — X, the second term in 
the brackets may be omitted: [X], [jl], and 01] stand for [X, VT - X], VT - X], 
and [!^, yP — X], respectively. Finally, for a function ip on the edges, we write y] 
(rather than ip{[X,Y])) for EeG[x,y] 

A tuple B = (r' I A, M; Bi, . . . , Bk), where T' is some BD-graph equivalent to F, is 
called an odd barrier for F if the following conditions hold with respect to F' (see Fig. [2]): 



(5.1) (i) A,M,Bi,...,Bk give a partition of VV = VT, and q e A. 

(ii) For each i = 1, . . . , k, c[^, Bi] is odd. 

(iii) For distinct i, j = 1, . . . ,k, c[Bi, Bj] = 0. 

(iv) For each i = l,...,k, c[Bi, M] = 0. 

The capacity of B is defined to be 

(5.2) c{B) ■.= 2c[t,% + c[l]-k. 



Theorem 5.1 (Max IBD-Flow^ Min Odd Barrier Theorem) ForT,c,q as above, 
the maximum IBD-flow value is equal to the minimum odd barrier capacity. A (feasible) 
IBD-flow g and an odd barrier B = (F' | A, M; Bi, . . . , Bk) for T have maximum value and 
minimum capacity, respectively, if and only if the following conditions hold with respect 
to T': 

(i) g[l, % = c[t, "a] and gfl, t] = 0; 
(h) g[l, M] = c[l, M] and g[A, M] = 0; 
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(iii) for each i = 1, . . . , k, either g[^, Bi] = c[^, -Bi] — 1 and g[A, Bi] = 0, or g[^, Bi] = 
c[l,Bi] andg[^,Bi] = 1. 

Note that there may exist many minimum capacity odd barriers for F, c, q. It is well- 
known that in a usual edge-capacitated digraph with a source s and a sink t, the set 
of nodes reachable by paths from s in the residual digraph of a maximum s-t flow / 
determines a minimum capacity s-t cut. Moreover, this minimum cut does not depend 
on the choice of / and, therefore, may be regarded as the canonical one. 

A similar phenomenon takes place for maximum IBD-flows and minimum odd barriers 
(and we will essentially use this in the proof of Proposition [33]) • To describe this, consider 
an IBD-flow g in T. The residual BD-graph Tg endowed with the residual capacity 
Cg : ETg is constructed in a similar way as for usual flows. More precisely, VTg = 

VT and the edges of Tg are as follows: 

(5.3) (i) each edge e G ET with g{e) < c(e) whose residual capacity is defined to be 
Cg{e) := c(e) - g{e), and 
(ii) the reverse edge to each edge e € ET with g{e) > 0; the directions of 
at the endpoints are reverse to those of e and the residual capacity is 
Cg{e^) ■■= g{e). 

A bidirected walk P in Tg is called Cg-simple if P passes each edge e at most Cg{e) 
times. If P is a Cg-simple closed q-q walk leaving its end q twice, we can increase the 
value of g in T by 2, by sending one unit of flow along P. So the existence (in Tg) of 
such a walk P, which is called {T,g) -residual, implies that g is not maximum. A converse 
property holds as well. 

Theorem 5.2 An IBD-flow g inT is maximum if and only if there is no (T, g) -residual 
walk. 

When we are given a maximum IBD-flow q, a certain minimum odd barrier can be 
constructed by considering the residual graph Tg. Namely, let Rr,g (resp. Rr,g) be the 
set of nodes v that are reachable by a (T, (7)-residual q-v walk that leaves q and enters v 
(resp. leaves both q and v). Then q ^ Rr,g, by the maximality of g. 

Theorem 5.3 Let g be a maximum IBD-flow for T, c, q. Define A := 
and M := VT - U R), where g and R := Rr,g- Let Bi, . . . , Bj^ be the node 

sets of weakly connected components of the underlying undirected subgraph ofTg induced 
by ^ n^. Define T' to be the BD-graph obtained from T by flipping the set li-t 
(contained in A). Then Bg := {T' \ A, M; Bi, . . . , B^) is a minimum odd barrier. 

An important fact is that the minimum odd barrier Bf does not depend on g, and we 
refer to it as the canonical odd barrier for T, c, q. 

Theorem 5.4 The sets ^r,g same for all maximum IBD-flows g in T, and similarly 
for the sets ^r,g, the minimum odd barriers Bf, and the graphs T' obtained from T by 
flipping ^r,g - '^r,g- 
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5.2 Proof of Proposition 13.31 



In fact, we have freedom of choosing any of the two forms (expensive or compact) of H 
to prove Proposition 13.31 in full, as it is easy to see that the claims in these cases are 
reduced to each other. We prefer to deal with the expensive form, taking advantage from 
certain nice structural features arising in this case. One reason for our choice is that any 
IBD-flow in the expensive H is automatically good, as explained in Remark 13.21 

We know that there exists a good fractional bidirected g-flow / in H that saturates 
the set £"0 of locked edges, and our goal is to show the existence of an IBD-flow saturating 
Eq. Recall that any edge e £ Eq is generated by some node v of G, i.e. e = e^. 

It will be convenient for us to construct the desired IBD-flow without explicitly im- 
posing the "lower capacities" on the locked edges. For this purpose, we modify H as 
follows. 

First, we add a loop V"^ with infinite capacity (entering q twice). Also we add to H 
a node z, which is regarded as a new source. 

Second, let Eq contain an edge = if'^lf''^ generated by a vertex v G VGi in some 
zone V^, s € T. We delete from H and, instead, add two edges if^V and ^if^ with 
capacity c{v) each. 

Third, let Eq contain the loops e^i {i = 1, . . . , c{w)) for some central node w oi Gi. 
We replace each e^i (having unit capacity) by edge z v~i with capacity 2; we call it the 
root edge at 9^i. 

We denote the resulting BD-graph by and keep the previous notation c for its 
edge capacities. The above g-flow / is transformed, in an obvious way, into a feasible 
z-flow in H^, denoted by / as before. Note that this / saturates all edges created from 
those in Eq (i.e. from e„ and e^i as above); these edges leave z and the value of / is 
maximum among the feasible z-flows in and is equal to 2c(-E'o)- 

Let g he & maximum IBD-flow in H^. We are going to prove that val(/) = val(g). 
This would imply that the corresponding IBD-flow in H saturates Eq as required. To 
this aim, consider the canonical odd barrier B = {H'^ \ A, M; Bi, . . . , Bj.) for i?^, c, z (see 
Theorem 15. 4p . Here is the BD-graph (with the source z) obtained from according 
to Theorem 15.31 (i.e. := F' for F := H"^)- From now on, speaking of edge directions, 
the capacities c and the flow g, we mean those in H^, unless explicitly stated otherwise. 

We have (cf. (lOD l 

(5.4) val(5) = c{B) = 2c[lt, ^4] + c[l] - k. 

The following assertion is crucial. 
Lemma 5.5 For each p = 1, . . . ,k: 

(i) Bp = {9j^,i} for some w EV''^ and z G {1, . . . , c{w)}; 

(ii) e^i is not locked (so contains the loop e^i but not the root edge at ); 

(iii) among the edges (legs) connecting A and Bp, one edge leaves A and the other edges 
enter A. 
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Proof. By the constructions of H and H^, for any w €zV^ and distinct i,j = 1, . . . ,c{w), 
there is an automorphism vr = '7Tw,i,j of that swaps O^^i and and is invariant on 
the other nodes. Also vr respects the capacities in H^, and the function g induced by g 
under tt (i.e. g{e) := g{7r{e))) is again a maximum IBD-flow in H^. Since ;B is canonical, 
it follows from Theorem 15.41 that 

(5.5) 'RHi,g=Rm,9, Rm,9 = ^m,g, Rm,9^Rm,g = AuBiU...uBk. 

The nodes in R^ig — '^^i g are flipped when constructing from . Then (|5.5p and 
Theorem 15.31 imply that for i,j as above, 

(5.6) (a) O^i is flipped if and only if O^j is flipped; 

(b) Oy^i G A if and only if O^j € A. 

(c) eBiVJ ...UBk'ii and only \i 6^^ G -Bi U . . . U Bk- 

Let p G {1, . . . , k}. Since the capacity c[X, Bp] (in H'^) is odd, the set [^, Sp] contains 
an edge e with c(e) odd. Any edge in having an odd capacity is either a loop e^i or 
a leg e^i g (regarding "infinite" capacities as even ones). 

Obviously, no loop can be "responsible" for the oddness of c[^. Bp]. 

So e = e^i g = 9^i9u,,s for some tt; G 1 < i < c{w) and s ^ T. Let e denote the 
edge of corresponding to e. Then (see Fig. [^b)) e leaves 6u},s and enters 9^i. Due to 
flips, however, this may not be the case for e in H^. 

Suppose G A (and 9u],s & Bp). Then e leaves 9^i, whence 9^i is a flipped 
node in A. Now (|5.6p (a.b) imply that all 9^j are flipped nodes belonging to A and that 

s ^ -^p] foi' all i = 1) • • • 1 c(ti;). But then e cannot be "responsible" for the oddness 
of c\~^, Bp] since c{w) is even. 

So we have 0^;^^ G ^ and G Bp. Then e leaves The edge e leaves ^^^^ as 

well. Hence 9uj^s is not flipped. Since c{w) is even, there must be j G {1, . . . , c{w)} such 
that the leg e^j ^ = 9w^s9^j is not in [^,i?p] (for otherwise one may pick another pair 
w',i'). Then 9^j is not in Bp. In view of (|5.6p (c). 9^j belongs to a S-set in B different 
from Bp. Considering the automorphisms vr = T^w,i',j' for all distinct = 1, . . . ,c{w) 
and using the fact that the canonical barrier B preserves under vr (in view of ()5.5p ). we 
can conclude that the nodes 9^i, . . . ,9^c(w) belong to different i?-sets in B. Since these 
S-sets are pairwise disjoint and each automorphism vr swaps two copies of 0^, and do 
not move the remaining nodes in H^, each of these i?-sets can contain only a single node. 
Thus, Bp = {9y^i}, yielding (i) in the lemma. 

Next we show (ii). From the construction of it follows that 

(5.7) A = ^H2g-%H2g and BiU . . .U Bk = ^H'^^g^^H'^,g- 

By the flrst equality, any (i/^, g)-residual walk ending at a node v (z A enters v, and by 
the second equality, there exist an (i?^, g')-residual walk P to 9 := 9^i that enters 9 and 
an (iJ^, g')-residual walk Q to 9 that leaves 9. Recall that the residual walks leave the 
source z. Let a = u 9 be the last edge of P, and b = v 9 the last edge of Q. Deflne E' 
{E") to be the set of legs e = e^i ^ with g{e) = (resp. g{e) = 1). Note that (cf. (|5.3p ) 
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(5.8) if e G £" then ^ EHg and e enters 9 in Hg, and if e G then e ^ E^F^ and 
leaves 6 in ff^. 

Supposing the existence of the root edge r = ^ 6 (in and i/^), we can come to a 
contradiction as follows. Since there is no loop at 6, both nodes u, v are in A. Note that 
the edge a is different from r (which leaves 6) and from (which enters z). Then ()5.8p 
implies that a £ E' . Furthermore, a is of the form it 6 . For if a enters u then the edge 
of P preceding a leaves u, whence the part of P from z to u forms an (i?^, (7)-residual 
walk leaving u, which is impossible since u (z A. 

So a G and g{a) = = c(a) — 1. Then g[A,Bp] = 0, by Theorem IS.lf iii). 

This implies that E" C [^,B]. But then the last edge b = v 9 of the walk Q as above 
cannot be reverse to any edge in E"; for otherwise b enters v, implying that the part of 
Q from z to V leaves v. Also b is neither reverse to an edge in E' (cf. (jS.Sp ). nor equal 
to r. The latter is because r G [3^,i?p], and therefore, g{a) < c{a) implies g{r) = c(r) 
(cf. Theorem I5.1( iii)). whence r ^ ^^g- Thus, Q does not exist. This contradiction 
yields (ii). 

It remains to show (iii). By (ii), we have e^i G EH"^ , and [A, Sp] is exactly the set 
of legs at 9 := 9^i. Suppose d := |[^,-Bp]| ^ 1. Then d > 3, since c[^,-Bp] = d is 
odd. Hence g[X,Bj^ > d—\ > 2 (by Theorem IS.ll fiii)). Also the fact that all legs 
enter 9 together with divg(0) = and gie^ji) < 1 implies that the only possible case is 
when g{e^i) = 1, g[^,Bp] = 2 and g[A,Bp] = 0. Now take an (ff^, (5r)-residual walk Q 
to 9 that leaves 9, and let b be its last edge. Then b is neither the loop e^i (which is 
saturated), nor reverse to a leg e = v 9 with g{e) > 0. Indeed, if b = then b enters v 
(in view of f G ^ and e G Bp]), and hence the part of Q from z to v leaves v, which 
is impossible since v G A. This contradiction yields (iii) and completes the proof of the 
lemma. □ 

Based on Lemma [5.51 we now finish the proof of Proposition |3]3l Consider Bp = {9^i} 
and let e^i., be the unique edge in Then [^4 , Sp] = {e^,^^^ \t£T—{s}}. Consider 

the maximum fractional BD-flow / as before. By the goodness of / (see ()3.2p ). we have 

f[l, Bp] - f^A, Bp] = f{e^.J - E,eT-{s} 

= /(e^,,J - (2/(e^.O - /(e^,,J) = 2 (/(e^.,J - /(e^.)) < = c[l,Bp] - 1. 
Using this and (15. 4p . we have 

val(/) = div^(z) = Y.^^^'^^Av) = (2f['t,% - 2f['A,l]) + (/[I] - fi^A]) 

= (2f[l, "A] - 2ffA, 1]) + (^f[l, M] - ffX, M]) + ^'^^ (/[t. Bp] - f[A, Bp]) 

< 2c[l, 'A]+c[1, + X]p_i (^[^' ^p] ~ = c[l] + 2c[l, 'A]-k = c{B) = V3l{g). 

Thus, we obtain the desired relation val(/) < val(g) (which, in fact, holds with equality). 
This completes the proof of Proposition 13.31 
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6 Dual half-integrality 



6.1 Polyhedral approach 

Theorem 6.1 Let a : VG —?■ Z+ and p G Z+. Then problem (D\) has a half-integer 
optimal solution. 

Proof. The proof follows easily from Theorem 12.11 and the general fact that the "to- 
tally dual l/A:-integrality" implies the "totally primal 1/ /c-integrality", which is a natural 
generalization of a well-known result on TDI-systems due to Edmonds and Giles |EG77| . 
More precisely, we utilize the following simple fact (see, e.g., fKa89, Statement 1.1]): 

Lemma 6.2 Let A be a nonnegative m x n-matrix, b an integral m-vector, and k a 
positive integer. Suppose that the program D{c) := max{y6 | y € Q!^, yA < c} has 
a 1/k-integer optimal solution for every nonnegative integral n-vector c such that D{c) 
has an optimal solution. Then for every nonnegative integral n-vector c, the program 
P{c) := minjcx | x G Q", Ax > b} has a 1/k-integer optimal solution whenever it has 
an optimal solution. 

In our case, we set k := 2 and take as A (resp. b) the constraint matrix (resp. the right 
hand side vector) of {D\). Then b is integral, D{c) becomes (A^a)) -P(c) becomes {D\), 
and the half-integrality for the former implies that for the latter. □ 

This proof is not "constructive" and does not lead directly to an efficient method for 
finding a half-integer optimal solution / to (-Da)- Below we devise a strongly polynomial 
algorithm. 

6.2 The algorithm 

It has as the input arbitrary (rational- valued) optimal solutions / and F to (Da) and (A'^a)) 
respectively, and outputs a half-integer optimal solution / to (-Da)- (Such / and F can be 
found in strongly polynomial time as described in Sectional) 

As before, we set i := a + 1, and in what follows, speaking of a geodesic, we always 
mean an ^-geodesic in G, i.e. a T-path P with i{P) = a{P) + 1{P) = A. Our goal is to 
construct I : VG satisfying the following conditions: 

(6.1) (i) a{P) +T{P) > A for any T-path P in G; 
(ii) a{P) + 1{P) = A for each geodesic P; 
(ih) for V G VG, if l{v) = then T{v) = 0. 

Then ()6.ip and the complementary slackness conditions ()2.6p - ()2.7p imply that the node 
length / forms an optimal solution to (-Da). 

We construct an undirected graph T and endow it with edge lengths n : ET as 
follows. We first include in T the terminal set T and all nodes and edges of G contained 
in geodesies. Also we add to T the edges of G with both ends lying on geodesies or T. 
The edges e of the current T are called regular and we define /i(e) := 0. 

Next we add to T additional edges, which are related to constraints due to parts of 
G outside F. More precisely, we scan all pairs of nodes u,v & VT not connected by a 
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(regular) edge and such that there exists a path Q in G having all nodes in VG — VT 
and whose first node is adjacent to u, and the last node to v. We add to F edge e = uv, 
referring to it as a virtual edge, and define its length ^{e) to be the minimum value of 
a{Q) among such paths Q. The construction of F, /x reduces to a polynomial number of 
usual shortest paths problems in G. 

Note that l{v) = holds for each node v eVG-VF (by and ^1}). We assign 
l{v) := for these nodes v and will further focus on finding values of I on the nodes in F. 

For a path P in F, let A(P) denote its full length a{P) + l{P)+fi{P). Clearly A(P) > A 
holds for any T-path P in F, and for each T-path Q in G, there exists a shortcut path P 
in F such that A(P) < a(Q) + 1{Q). 

The desired lengths / on VT will be extracted from a system of linear constraints 
described below. For a node v € ^F, let Ty (Jly) denote the set of terminals s € T (resp. 
pairs s,t £ T) such that v belongs to a geodesic from s (resp. connecting s and t). When 
a terminal s belongs to no geodesic, we set by definition := {s}. For each v € VT and 
s G Ty, we introduce two variables pj{v) and pf{v) and impose the following constraints: 



(6.2) (i) For each s e T, p7(s) = 0. 

(ii) For each v € VT and s (zTy, 

pt{v)-p-{v) =a{v) if l{v) = 0, 
> aiv) if l{v) > 0. 

(iii) FoT each V G VT and {s,t} £ Tly, pf{v) + p^{v) = X{andpf{v)+pj{v) = X). 

(iv) li e = uv £ ET and s G TyDTy, then 

Pi{v)-Ptiu) <p{e), 
Pi{u)-pt{v) </i(e). 

Moreover, if there exists a geodesic from s containing both u, v in this order 
(resp. in the order v,u), then the former (resp. the latter) inequality is 
replaced by equality. (Note that in this case p{e) = 0.) 

(v) If e = uf G ET, s E Ty, t G Ty, and s ^t, then pf{u) + pf{v) > A — p{e). 
The meaning of these variables becomes evident from the proof of the next statement. 



Lemma 6.3 System i6.2\) has a solution. 

Proof. For a p-q path P in F, define its pre-length to be A(P) — (a(g) + l{q)) (i.e. 
compared with the full length, we do not count the last node). For v € VT and s S Ty, 
define pj{v) (resp. pf(v)) to be the minimum pre-length (resp. the minimum full length) 
of an s~v path in F. Then (|6.2p (i)-(iii) follow from the construction. Condition (|6.2p (iv) 
represents a sort of triangle inequalities (giving one equality if e belongs to a geodesic 
from s). Finally, condition (|6.2p (v) holds since the full length of any s-t path in F is at 
least A. □ 
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We observe that in linear system ()6.2p . each constraint contains at most two variables, 
each occurring with the coefficient 1 or -1, and that the R.H.S. in it is an integer. A 
well-known fact is that a linear system with such features is totally dual half-integral; 
therefore, it has a half-integer basis solution (whenever it has a solution at all), and such 
a solution can be found in strongly polynomial time (cf., e.g., |EJ70| [5c03| ). 

Given a half-integer solution (p^ , p^) to ()6.2p . we define half-integer node lengths / 
as follows: 

l{v) := p1i'{v) — pj{v) — a{v) for all v € VT and s € T^. 
Now the desired algorithmic result is provided by the following 

Lemma 6.4 / is well-defined and satisfies \6.1\) . 

Proof. We first show that for any v G VT and s,t G T^, 

(6-3) pt{v) - Pi{v) = pt{v) - pt{v). 

This is trivial when = (since in this case v G T and T„ = {v}). Let 11^, ^ 0. 
If {s,t} G n^, then ([OD follows from ([H2])(iii)- Now with any two s,t e % is 

implied by the fact that the graph whose nodes and edges are the elements of and Yiy , 
respectively, is connected (as it is easy to see that for g} G H^, at least one of 

{s,p}, {s,q\ is in 11^, as well). So / is well-defined. 
Property (|6.1|) fiii) is immediate from (|6.2|) fii). 

To see ()6.ip fii). consider an s-t geodesic P. Going along P step by step and apply- 
ing (|6.2p (ii).(iv). we observe that for each node v on P, the s~v part P' of P satisfies 
a(P') +1{P') = pt{v)-p~{s). When reaching t, we obtain a{P) +T{P) = p+{t)-pj{s), 
and now follows from ([621) (iii) and pj{s) = p~[ {t) = (by ([621) (i))- 

Finally, consider an arbitrary T-path Q in T, from p to q say. To conclude with (|6.ip (i). 
it suffices to show that 

(6.4) A(Q) := a{Q) + T{Q) + p{Q) > A. 

Represent Q as the concatenation Q' ■ Q" , where Q' is a part of a geodesic from p. 
We prove ()6.4p by induction on the number \Q"\ of edges in Q" . When \Q"\ =0, Q is 
a geodesic, and we are done. Assuming this is not the case, take the first edge e = uv 
of Q" , where u is the end of Q' . By reasonings above, A(Q') = Pp{u). If u G T (and 
therefore, v = q), ()6.4p immediately follows from ()6.2p (v) (with s := q and t ■= q)- And 
a V ^ T, then v belongs to some s-t geodesic L. W.l.o.g., one may assume that s p 
and t ^ q. Applying ()6.2p (v) to s,p, e, we have 

P^{u) + pt{v)+p{e)>X. 

Comparing this with pj{v) + Pt{v) = A and using p~^{v) — pj{v) = a{v) + l{v), one can 
conclude that A((5) > A(i?), where R is the t-q path being the concatenation of the t~v 
part of (the reverse of) L and the v-q part R" of Q. Since = \Q"\ — 1, we can apply 
induction and obtain A{Q) > A{R) > A, as required. □ 
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Appendix: Skew-symmetric graphs and flows 



In this section we recall the notions of skew-symmetric graphs and integer skew-symmetric 
flows, review known results on such graphs and flows, and then use them to derive 
necessary results on bidirected graphs and flows to which we appealed in Section [5l 

7.1 Skew-symmetric graphs 

A skew- symmetric graph, or an SK-graph for short, is a digraph G = {VG, EG), with 
possible parallel arcs, endowed with two bijections av,(JA such that: ay is an involution 
on the nodes (i.e. av{v) ^ v and av{crv{v)) = v for each node v); a a is an involution 
on the arcs; and for each arc a from u to v, 0"^(a) is an arc from (yy{v) to ay(u). 
For relevant results on SK-graphs and a relationship between SK- and BD-graphs, see 
[Tu67| IGK96| IGK04| IBK07| . For brevity ay and a a are combined into one mapping 
a on VG U AG, which is called the symmetry (or skew-symmetry, to be precise) of G. 
For a node (arc) x, its symmetric node (arc) 0"(x) is also called the mate of x, and we 
usually use notation with primes for mates, denoting ci(x) by x' . Although G is allowed 
to contain parallel arcs, when it is not confusing, an arc from uto v may be denoted as 
(n, v) or u^. 

Observe that if G contains an arc a from a node v to its mate v' , then o! is also an 
arc from v to v' (i.e. a' is parallel to a). 

The symmetry a is extended in a natural way to walks, subgraphs and other objects 
in G. In particular, two walks are symmetric to each other if the elements of one of 
them are symmetric to those of the other and go in the reverse order: for a walk P = 
{vo,ai,vi, . . . ,ak,Vk), the symmetric walk P' = a{P) is (w^, a'f^,v'^_^, . . . ,ai,Uo). 

Next we explain a relationship between skew-symmetric and bidirected graphs. Given 
an SK-graph G, choose an arbitrary partition vr = {Vi, V2) of VG such that V2 = c(Vi). 
Then G and vr determine the BD-graph G* with VG* = Vi whose edges correspond to 
the pairs of symmetric arcs in G. More precisely, arc mates a, a' of G generate one edge 
e of G* connecting nodes u,v ^Vi such that: (i) e goes from u to v \i one of a, a' goes 
from n to v (and the other goes from v' to u' in V2)', (ii) e leaves both u, v if one of a, a' 
goes from u to v' (and the other from v to u'); (iii) e enters both u, v if one of a, a' goes 
from u' to V (and the other from v' to u). Note that e becomes a loop if a, a' connect a 
pair of symmetric nodes. 

Conversely, a BD-graph G* determines an SK-graph G with symmetry a as follows. 
Make a copy (t{v) of each element v of V* := VG*, forming the set (V*)' := {(t{v) \ 
V G V*}. Put VG := V* U {V*y. For each edge e of G* connecting nodes u and v, 
assign two "symmetric" arcs a, a' in G so as to satisfy (i)-(iii) above (where u' = cr{u) 
and v' = cr{v)). An example is depicted in Fig. [3j 

Remark 7.1 Note that one BD-graph generates one SK-graph, by the second construc- 
tion. On the other hand, one SK-graph generates a set of BD-graphs, depending on the 
partition vr of V, by the first construction. Namely, for each pair of symmetric mates 
{v,v'} in G one may distribute v,v' between Vi,V2 so that either v G Vi, € V2 or, re- 
versely, V G V2, v' GVi. The resulting BD-graphs are obtained from one other by making 
corresponding flips (defined in Subsection \5. 
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(a) BD-graph G*. (b) Corresponding SK-graph G. 

Figure 3: Related bidirected and skew-symmetric graphs. 



There is essentially a one-to-one correspondence between the walks in G* and G. 
More precisely, let r be the natural mapping of VG U AG to VG* U EG* . Each walk 
P = {vq, ai,vi, . . . , ak,Vk) in G (where Oj = (uj-i, Wj)) induces the sequence 

t{P) ■= {T{vQ),T{ai),T{vi), . . . ,T{ak),T{vk)) 

of nodes and edges in G* . One can see that t[P) is a walk in G* (i.e. T(ai), T(aj_|_i) form 
a transit pair at T{vi), for each i) and that t{P') is the walk reverse to t{P). Moreover, 
for any walk P* in G* , there is exactly one walk P in G such that t{P) = P* (considering 
P up to replacing an arc a € AP by its mate a' when a, a' are parallel, i.e. correspond 
to a loop in G*). 

7.2 Skew-symmetric flows 

We call a function ip on the arcs of an SK-graph G (self-) symmetric if '^{a) = '^{a') for 
all a G AG. Let s G VG be a designated source; its mate s' is regarded as the sink. 
An integer skew-symmetric s-s' flow, or an ISK-flow for short, is a symmetric function 
/: AG — > Z+ being an s-s' flow in a usual sense: div f{v) = for all v E T^G — {5,5'}, 
and divj(s) > 0. The value of / is val(/) := divj(s). Here <l\Vf{v) denotes the usual 
divergence (given by (|3.1|) . where 5™(?;) and 5°^^{v) are the sets of arcs entering and 
leaving u, respectively). 

For a capacity function c: AG — >■ Z+, a flow / is said to be feasible if /(a) < c(a) for 
all a € ^G. We refer to a feasible ISK-flow of maximum possible value as a maximum 
ISK-flow. 

The above correspondence between BD- and SK-graphs is naturally extended to flows. 
More precisely, if / is a symmetric ss' flow in G, then transferring the values of / from 
the pairs of arc mates of G to the edges of the BD-graph G* := t(G), we obtain a 
r(s)-flow in G*, denoted as /*. The converse correspondence is evident as well. 

For X, y C VG, let {X, Y) denote the set of arcs going from X to Y . Also (ac- 
commodating notation from Section [5] to digraphs) we denote by \^] the set of arcs 
leaving X. 
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Figure 4: A skew-symmetric odd-barrier. Grayed arcs correspond to odd capacity con- 
straints. 



Let c : AG — t- Z^. be a symmetric capacity function. Then a tuple ,B = 
(A, M; . . . , of subsets of is called a (skew-symmetric) odd harrier (w.r.t. 
the source s) if the following conditions hold (see Fig. 

(7.1) (i) the sets A^A! = a{A), M, Bi, . . . ,Bk give a partition of VG, each Bi is self- 

symmetric (B'- = Bi), and s A; 

(ii) For each i, c{A,Bi) is odd. 

(iii) For distinct c{Bi,Bj)=0. 

(iv) For each i, c{Bi,M) = c{M,Bi) = 0. 

The capacity of B is defined to be 

(7.2) c(;B) := cil] - k. 

Odd barriers in skew-symmetric graphs are related to their bidirected counterparts 
introduced in Section O Indeed, consider a BD-graph G* with integer edge capacities 
c: EG* — > Z+ and a source s. Construct the related SK- graph G with VG = V UV' , 
where V := VG*. Edge capacities c in G* induce symmetric arc capacities in G, also 
denoted by c. The source s in G* gives the source s and the sink s' in G. Consider a 
skew-symmetric odd barrier B = {A, M; Bi, . . . , B^) in G. 

This barrier gives rise to the following odd BD-barrier B* in G* obeying c{B*) = c{B). 
We first construct a new BD-graph from G by taking a bipartition (Vi, V2 = aiVi)) of 
yG such that ACVi and - (^ U A') = - U A'); cf. Remark O The resulting 
BD-graph H* is equivalent to G*. Moreover, i?* is obtained from G* by flipping a subset 
of nodes within A. 

The node subsets M, Bi, . . . , Bf^ in G are self-symmetric and induce subsets 
M*,Bl,. . . ,Bl in G* and H* in a natural way; namely, M* := M n F = M n 
and similarly for B* . Define ^* := {H*\A* ,M*;Bl, S*), where A* := (A U yl') n V. 
A straightforward examination shows that the properties in (|7.1|) imply their bidirected 
counterparts in ([SJ]). To see that c{B*) = c{B), define Z := M U Bi U . . . U Bk- Note 
that c[2] = c{A, A') + c{A, Z). The capacity c{A, A') is equal to 2c[A^, ^] (in H*) since 
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{A, A') consists of pairs of arc mates, each pair corresponding to an edge in [A , A*]. 
And the capacity c{A, Z) is equal to c[A^] since the (symmetric) set Z corresponds to 
M*\JBl...Bl in H*. 

In hght of these observations, Theorem 15.11 is a consequence of the fohowing Tutte's 
theorem. (For shorter proofs of this and next theorems, see also |GK04) .) 

Theorem 7.2 (Max ISK-Flow Min Odd Barrier Theorem |Tu67| ) ForG,c,s as above, 
the maximum ISK-flow value is equal to the minimum odd barrier capacity. An ISK- 
flow f and an odd barrier B = {A, M; Bi, . . . , Bk) have maximum value and minimum 
capacity, respectively, if and only if the following hold: 

(i) f{A, A'UM) = c{A, A' U M) and f{A' UM,A) = 0; 

(ii) for each i = 1,... ,k, either f {A, Bi) = c{A,Bi)-l andf{Bi,A) = 0, or f{A,Bi) = 
c{A,Bi) andf{Bi,A) = 1. 

Next we establish additional correspondences. Consider an ISK-flow / in G. The 
residual SK-graph Gf endowed with the residual capacities cj : AG — >■ Z_|_ is constructed 
in a standard fashion: VGf = VG, and the arcs of Gj are: 

(7.3) (i) each arc a G AG with f{a) < c{a) whose residual capacity is defined to be 
Cf{a) := c(a) — /(a), and 
(ii) the reverse arc = {v,u) to each arc a = {u,v) € AG with /(a) > 0; its 
residual capacity is Cf{a^) := /(a) 

(cf. (|5.3|) ). A path P in Gj is called Cf -regular if cj(a) = Cf{a') > 2 holds for each pair 
of arc mates a, a' occurring in P. (In other words, the bidirected image of P in is a 
cj-simple walk.) If P is a cj-regular s-s' path, we can increase the value of / by 2 (by 
sending one unit of flow along P and one unit of flow along P'). So the existence of such 
a P implies the non-maximality of /. A converse property is valid as well. 

Theorem 7.3 ( [Tu67 j) An ISK-flow f is maximum if and only if there is no Cf -regular 
s-s' path in Gf. 

This implies Theorem 15.21 for IBD-flows. 

Given a maximum ISK-flow /, a certain minimum odd barrier can be constructed 
by considering the residual graph Gj. The construction described in the proof of Theo- 
rem 3.5 in |GK04| (relying on Lemma 2.2 in [ GK96| ) is as follows. 

Theorem 7.4 Let f be a maximum ISK-flow. Let R = Rf be the set of nodes reachable 
from s by Cf-regular paths in Gf. Define A := R — R' and M := VG — {RU R'). Let 
Bi, . . . , Bk be the node sets of weakly connected components of the subgraph G induced by 
R n R' . Then Bf := {A, M; Bi, . . . , B^) is a minimum odd barrier. □ 

This subset R of nodes in G corresponds to two sets = '^G*,f* a-i^d ^ = Rq*j* in 
G* (defined just before Theorem 15. 3t here T = G* and g = /*). More precisely, assuming 
that each node v G VG* corresponds to node mates v,v' in G (cf. Subsection 17. ip . one 
can realize that (resp. R) is the set of nodes v € VG* such that v G R (resp. v' £ R). 

Finally, the last theorem in Subsection 15.11 (Theorem 15. 4p is implied by the following 
assertion. 
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Theorem 7.5 The sets Rf in Theorem \ 7.4\ are equal for all maximum ISK-flows f. 
Therefore, the minimum odd barriers Bf are equal as well. 

This fact can be extracted from reasonings in | iGK04) . yet it is not formulated there 
exphcitly. For this reason, we give a direct proof. 

Let / be a maximum ISK-flow such that the set Rf is inclusion- wise minimal and let 
Bf := {A, M; Bi, . . . , Bi;.). Consider another maximum ISK-flow g (if any). Then 

(7.4) c[t] - k = c{B) = val(5) = Y.,^^ div,(^;) 

= g{A, A') + g{A, M) + g{A, Bi) + . . . + g{A, B^), 

where for disjoint subsets X,Y C VG, g{X,Y) denotes g(X,Y) — g{Y,X). For i = 
1, . . . ,k, we have: g{A, Bi) < c{A, Bi); c{A, Bi) is odd; and g{A^ Bi) is even (the latter 
is due to a result in |Tu67| : see also |GK04| Corollary 3.2]). Therefore, c{A,Bi) — 
g{A,Bi) > 1. Also g{A,A') < c{A,A') and g{A,M) < c{A,M). Comparing these 
relations with (|7.4|) . we conclude that: 

(i) all arcs in {A, A' U M) are saturated by g, while all arcs a in (A' U M, A) are free 
of g (i.e. g{a) = 0); 

(ii) for each i, g{A,Bi) = c{A,Bi) — 1. 

In terms of the residual graph Gg, (i) and (ii) mean that the sets A and VG — A 
are connected in Gg by exactly k arcs ai , . . . , , each Oj going from A to Bi and having 
the residual capacity 1. By symmetry, A' and VG — A' are connected by only the arcs 
a'l, . . . ,a'f^ (each goes from Bi to A', and Cg{ai) = 1). Also by (|7.1|) fiii).(iv). no arc in 
Gg connects different sets among M,Bi, . . . ,Bk. Therefore, the set Rg of nodes in Gg 
reachable from s by c^-regular paths is contained in Rf. By the minimality Rf, we 
have Rf = Rg, as required. □ 

Acknowledgements. We thank the anonymous referee for useful suggestions. 
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